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Abstract 

The appearance of infinitely-many period-doubling cascades is one 
of the most prominent features observed in the study of maps de- 
pending on a parameter. They are associated with chaotic behavior, 
since bifurcation diagrams of a map with a parameter often reveal a 
complicated intermingling of period-doubling cascades and chaos. 

Period doubling can be studied at three levels of complexity. The 
first is an individual period-doubling bifurcation. The second is an 
infinite collection of period doublings that are connected together by 
periodic orbits in a pattern called a cascade. It was first described by 
Myrberg and later in more detail by Feigenbaum. The third involves 
infinitely many cascades and a parameter value (jq of the map at which 
there is chaos. We show that often virtually all (i.e., all but finitely 
many) "regular" periodic orbits at ^2 are each connected to exactly 
one cascade by a path of regular periodic orbits; and virtually all cas- 
cades are either paired - connected to exactly one other cascade, or 
solitary - connected to exactly one regular periodic orbit at f/,2- The 
solitary cascades are robust to large perturbations. Hence the inves- 
tigation of infinitely many cascades is essentially reduced to studying 
the regular periodic orbits of F(fi2, •)• Examples discussed include the 
forced-damped pendulum and the double-well Duffing equation. 

1 Introduction 

In Figure [TJ as fi increases from fi = —0.25 towards a value \if ~ 1.4, a fam- 
ily of periodic orbits undergoes an infinite sequence of period doublings with 
the periods of these orbits tending to 00. This infinite process is called a cas- 
cade. We will later define it more precisely. It has been repeatedly observed 
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Figure 1: Cascades for F(fi,x) = fi — x 2 . This figure shows the attracting 
set for F for —0.25 < /i < 2. The attracting set is created at a saddle-node 
bifurcation at \i = —0.25 (green dot). The path of unstable fixed points 
(red) exists for all fj, > —0.25. The stable fixed point undergoes infinitely 
many period-doubling bifurcations, limiting to the value \i ~ 1.4. This set 
of period doublings is called a period-doubling cascade. This map also has 
infinitely many period-doubling cascades that begin with periodic orbits of 
period > 1. The red curve consists of unstable regular fixed points that exist 
for all fi > -0.25. 



in a large variety of scientific contexts that the presence of infinitely many 
period-doubling cascades is a precursor to the onset of chaos. For example, 
cascades occur in what numerically appears as the onset to chaos for both the 
double- well Duffing equation, as shown Figure [2j and the forced-damped pen- 
dulum, shown Figure [3] Cascades were first reported by Myrberg in 1962 [9], 
and studied in more detail by Feigenbaum [4J. This cascade is not the only 
cascade. In fact, there are infinitely many distinct period-doubling cascades. 
Namely, there are infinitely many windows, that is, disjoint intervals in the 
parameter that begin with a saddle-node (or source-sink) bifurcation, and 
continue with the attractor undergoing an infinite sequence of period dou- 
blings within that interval of parameters. 
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Figure 2: Cascades in the double-well Duffing equation. The attract- 
ing sets (in black) and periodic orbits up to period ten (in red) for the time-27r 
map of the double-well Duffing equation: x"{t) + 0.3x'(£) — x{t) + (x(t)) 2 + 
(x(t)) 3 = /isint. Numerical studies show regions of chaos interspersed with 
regions without chaos, as in the Off-On-Off Chaos Theorem (Theorem [5]). 




Figure 3: The forced-damped pendulum. For this figure, periodic points 
with periods up to ten were plotted in red for the time-27r map of the forced- 
damped pendulum equation: x"(t) + 0.2x'(t) + sm(x(t)) = /icos(t), indicating 
the general areas with chaotic dynamics for this map. Then the attracting 
sets were plotted in black, hiding some periodic points. Parameter ranges 
with and without chaos are interspersed. 
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Quadratic maps as in the example in Figure [T] has the quite atypical 
property that as the parameter increases, there are no bifurcations that de- 
stroy periodic orbits. Such maps are called monotonic. (This monotonicity 
was originally proved implicitly by Douady and Hubbard in the complex an- 
alytic setting. See the Milnor-Thurston paper [7j for a proof.) Once one 
knows that a map is monotonic, it is easy to show that as chaos develops 
there must be infinitely many cascades. See Figures [TJ and |4]-[8j 

The monotonicity property is a quite severe restriction, even in one- 
dimension. No higher dimensional maps that develop chaos are monotonic; 
yet numerical studies indicate that there are infinitely many cascades when- 
ever there is one. See for example Figures |2j [3| and [5] In this paper we 
summarize our progress and give new extensions to our theory that explains 
why there are infinitely many cascades in the onset to chaos. Our explanation 
is also valid for maps of arbitrary dimension. 

In the first result of this paper, we consider the context in which virtually 
all periodic saddles have the same unstable dimension. (By virtually we 
mean all except for a finite number.) In this case, the onset of chaotic 
behavior always includes infinitely many cascades. 

There is an extensive literature on Routes to Chaos; that is, situations 
in which for some /ii and /i 2 , the trajectory of x under F(fj, 2 , •) is in the 
basin of a chaotic attractor, whereas under F(fii, •) it is not. Whatever 
might have happened that caused this change between n\ and /i 2 is called a 
route to chaos. We prefer to call these "routes to a chaotic attractor" to be 
more specific. There are many different routes to a chaotic attractor. See 
our discussion section for a partial enumeration. For many maps there is 
competition between instability and stability. For example, the appearance 
of an attracting periodic orbit as a parameter is varied may mask the chaotic 
dynamics, and when the orbit becomes unstable, a chaotic attractor is likely 
to appear. Hence a periodic orbit's loss of stability is one example of a route 
to a chaotic attractor. This approach ignores the question we address here: 
how did the chaotic dynamics arise in the first place? 

Two types of cascades. For maps with the monotonicity property, 
each cascade is solitary, in that it is not connected to another cascade by a 
path of regular periodic orbits. These paths are the colored stems shown 
in Figure |4j See Section [2] for full definitions of these terms. Furthermore, 
the chaos persists for all parameters larger than a certain value. However, 
in many scientific contexts, it is common to see chaotic behavior appear and 
then disappear as the parameter \i increases. Thus as it increases there is 
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Figure 4: Cascades for F(fi, x) = fix(l — x). The logistic map has infinitely 
many cascades of attracting periodic orbits, and all cascades start at the 
stable orbit of a saddle-node bifurcation. The unstable orbits form what we 
call the stems of the cascades (shown in color). Each stem continues to exist 
for all large fi values. By our terminology, this means that all the cascades 
shown are solitary (on any parameter interval [/Xi,/i 2 ], for /i X = 3.5 and any 
^2 > 4) since the stem does not connect its cascade to a second cascade. The 
stems are shown here up to period six. Different colors are used for different 
periods. 



both a route to chaos followed by a route away from chaos. In this situation 
virtually all cascades are paired; that is, two cascades are connected by a 
path of regular periodic orbits. See Figure |5j 

Solitary cascades are robust. In our second set of results, we show 
that while paired cascades can be easily created and destroyed, solitary cas- 
cades are far more robust even in the presence of rather large perturbations 
of a map. Solitary cascades usually have stems with a constant period. This 
stem-period can be thought of as the period that starts the cascade. These 
ideas give quite striking results. For example for each period p the following 
two maps 

Q(fi, x) = fi — x 2 

and 

Q(fx, x) = jj, — x 2 + 1000 cos(/i 3 + x) 
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Figure 5: Paired cascades in the Henon map (u, v) t— > (1.25— u 2 +fiv , u). 
The top bifurcation diagram shows a set of four period-7 cascades. The 
bottom bifurcation diagram shows detail of the top part. Only one point 
of each of the period-7 orbits of the Henon map are shown so that it is 
clearer how the two pairs connect to each other. The leftmost and rightmost 
cascades form a pair that is connected by a path of unstable regular periodic 
orbits (shown in red). Likewise, the two middle cascades form a pair. It is 
connected by a path of attracting period-seven orbits (blue). Paired cascades 
are not robust to moderate changes in the map. 
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have exactly the same number of solitary cascades of stem-period p - as- 
suming the bifurcations of the second map are generic. Namely, we call a 
map generic if its periodic orbit bifurcations are all generic. See Section [2] for 
the full definitions of these terms. We know that all the bifurcation orbits 
are generic for almost every smooth map, and that if a smooth map is not 
generic, then it has infinitely many generic maps arbitrarily close to it, but 
unfortunately - with few exceptions - we cannot tell if a given map is generic. 

The map Q has no paired cascades, but the Q may. For example there is 
exactly one solitary cascade with stem-period 1 and one with stem-period 3. 
These results extend to 

F(fj,, x) = fi — x 2 + g(fi, x) 

where g is smooth (ie. infinitely continuously differentiable) and \g(p, x) \ and 
\g x (n, x)\ are uniformly bounded - as would be the case if g was smooth and 
periodic in each variable, again assuming the map F has generic bifurcations. 

Outline. The paper proceeds as follows: In Section [2j we give some basic 
definitions, including what we mean by a cascade, the definition of chaos that 
we use here, and the class of generic maps with which we work. Section [3] 
contains a series of results relating chaos and cascades, with an explanation 
of the concrete relationship between periodic orbits within the chaos and 
the resulting cascades along the route towards this chaos. In Section [4], we 
discuss the fact that all cascades are either paired or solitary, and show that 
solitary cascades are robust under changes in the map. 

In Section [5j we show that if there is chaotic behavior interspersed with 
non-chaotic behavior, then virtually all cascades are paired. It is common 
in scientific applications that chaos is interspersed with orderly behavior, in 
what we call off-on-off chaos (defined formally in Section [5]). Our numerical 
studies indicate that this occurs multiple times for both the forced-damped 
pendulum and double-well Duffing examples. 

We end with a discussion and present open questions in Section [6j 



2 Definitions 

We investigate smooth maps F(fi, x) where \x is in an interval J, and x is in 
a smooth manifold SOT of any finite dimension. For example, for the forced 
damped pendulum, 

-j-t + A— + sin d = fi cos t, 
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Figure 6: Center manifold of saddle-node bifurcations and period- 
doubling bifurcations. This figure shows typical saddle-node and period- 
doubling (halving) bifurcations, along with the stability. For generic maps 
F : R x R n , it is sufficient to examine the R x R center manifold. We plot 
one-dimensional x vertically and \i horizontally. We use vertical arrows to 
show how the stability varies near a bifurcation periodic point - indicated 
by a large dot. In each case, all the stability arrows can be reversed, thereby 
generating four more cases. 



we take F(fi, x) to be the time-T map, where T = 27r is the period of the 
forcing, and x = (6, dO/di). Then the first coordinate of i is on a circle, and 
the second is a real number. Hence 9JT is a cylinder. 

We say a point (/i, x$) is a period-p point if F p (fi, xo) = xo and p is the 
smallest positive integer for which that is true. Its orbit, sometimes written 
[(//, Xo)], is the set 

{O,x ),(/i,xi),--- where xj = F 3 (//, x ). 

By the eigenvalues of a period-p point ((j,,xq), we mean the eigenvalues 
of the Jacobian matrix D x F p (^i,xq). 

An orbit is called hyperbolic if none of its eigenvalues has absolute value 
1. All other orbits are bifurcation orbits. Figure [6] depicts two standard 
examples of bifurcation orbits and the resulting stability of nearby periodic 
points. 

We call a periodic orbit a flip orbit if the orbit has an odd number of 
eigenvalues less than -1, and -1 is not an eigenvalue. (In one dimension, 
this condition is: derivative with respect to x is < — 1. In dimension two, 
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Figure 7: The regular periodic orbits form a one-manifold near 
regular saddle-node and period-doubling bifurcation orbits. In this 
schematic figure each point is an orbit and the horizontal axis is the parame- 
ter, usually /i in this paper, (a & b): Near a standard saddle-node bifurcation 
of a periodic orbit, the local invariant set consists of a curve of periodic orbits. 
They are either all flip orbits or all regular periodic orbits. Therefore RPO is 
locally a curve, (c & d): Near a standard period-doubling (or period-halving) 
bifurcation of a periodic orbit, the local invariant set consists of two curves 
of periodic orbits, one of period p shown as segment AB, and one of period 
2p shown as segment C. Segment C always consists of regular periodic orbits, 
whereas exactly one of A and B consists of flip orbits and the other regular 
periodic orbits. Thus RPO is locally a curve consisting of C and either A 
or B, depending on which is regular. For the quadratic map \x — x 2 , only 
(a) and (d) occur. That is, periodic orbits are created but never destroyed 
as \x increases. When x is two-dimensional, such simplicity virtually never 
occurs [5]. 



flip orbits are those with exactly one eigenvalue < — 1.) All other periodic 
orbits are called regular. For example, the periodic orbits of constant period 
switch between flip and regular orbits at a period-doubling bifurcation orbit 
since an eigenvalue crosses —1. See Figure [7| We write RPO for the set of 
regular periodic orbits. 

For some a,b G R and ip £ [ a >b), let Y(ip) = (fi(i[)),x(if))) be a path of 
regular periodic points depending continuously on ip. Assume ip does not 
retrace any orbits. That is, each Y(ip) is a periodic point on a regular periodic 
orbit, and distinct values of ip correspond to periodic points on distinct orbits. 
We call a regular path Y{ip) a cascade if the path contains infinitely many 
period-doubling bifurcations, and for some period p, the periods of the points 

in the path are precisely p, 2p, 4p, 8p, As one traverses the cascade, the 

periods need not increase monotonically, but as ip — > b, the period of Y{jp) 
goes to oo. The orbits of a cascade with monotonic period increase are 
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Figure 8: A depiction of a monotonic cascade. A cascade is a path of 
regular periodic orbits that has infinitely many period-doubling bifurcations 
with the periods going to infinity at the end of the path. Solid lines denote 
regular orbits, and dashed lines denote flip orbits. This figure uses only the 
orbit-creation bifurcations, (a) and (d) in Figure [7j If only bifurcation types 
(a) and (d) are present and x is scalar, then each saddle-node spawns an 
unstable branch and an attractor as /i increases. For fi — x 2 , there are no 
attractors for fx > // 2 for sufficiently large. Hence this branch of attractors 
cannot continue forever and must bifurcate. Only (d) is available, spawning 
a flip orbit branch and a doubled-period branch of attractors. Again this 
new branch of attractors cannot continue forever as \x increases. In this way, 
an infinite set of period-doublings results before /i reaches \xi. Much more 
complicated patterns are possible when all four types of bifurcations are 
allowed, including period-halving as well as period-doubling, and the path 
may contain many saddle-node bifurcations. 
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depicted schematically in Figure |8j 

Write fixed(fx,p) for the set of fixed points of F p (/i, •) and \fixed(/i,p)\ 
for the number of those fixed points. We say that there is exponential 
periodic orbit growth at fi if there is a number G > 1 for which the 
number of periodic orbits of period p satisfies \fixed(p,,p)\ > G p for infinitely 
many p. For example, this inequality might hold for all even p, but for odd 
p there might be no periodic orbits. This is equivalent to the statement that 
for some h(= logG) > 0, we have 

log|/ia;ed(/x,p)| 
hmsupp^ >h. (1) 

Periodic orbit chaos. Chaotic behavior is a real world phenomenon, and 
trying to give it a single definition is like trying to define what a horse is. 
Definitions are imperfect. A child's definition of a horse might be clear 
but would be unsatisfactory for a geneticist (whose definition might be in 
terms of DNA) and neither would satisfy a breeder of horses who might give 
a recursive definition, "an offspring of two horses." Definitions of real-life 
phenomena describe aspects of that phenomena. They might agree in the 
great majority of cases on which animals are horses, though there may be 
rare atypical exceptions like clones where they might disagree. Just as it 
is impossible currently to connect the shape and sound of a horse with its 
DNA sequence, it is similarly impossible currently to identify in full generality 
positive Lyapunov exponents with exponential growth of periodic orbits. 

Similarly "chaos" and "chaotic" should have definitions appropriate to 
the needs of the user. On the other hand, an experimenter might insist that 
to be chaotic, there must be a chaotic attractor, until he/she starts looking 
for chaos on basin boundaries and finds transient chaos. That approach 
leaves no terms for the chaos that occurs outside an attractor, as on fractal 
basin boundaries. We make no such restriction. Our results involve periodic 
orbits, and we make our definition accordingly. 

We say a map F(fi, •) has periodic orbit (PO) chaos at a parameter /i 
if there is exponential periodic orbit growth. This occurs whenever there is 
a horseshoe for some iterate of the map. It is sufficiently general to include 
having one or multiple co-existing chaotic attractors, as well as the case 
of transient chaos. As hinted at by Equation [TJ in many cases PO chaos 
is equivalent to positive topological entropy. We discuss this relationship 
further in Section [U 
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The unstable dimension Dim u (//, Xq) of a periodic point (//, x ) or pe- 
riodic orbit is defined to be the number of its eigenvalues A having |A| > 1, 
counting multiplicities. 

We say there is virtually uniform (PO) chaos at fi if there is PO 
chaos, and all but a finite number of periodic orbits have the same unstable 
dimension, denoted Dim u (/x). 

For the pendulum map discussed above, whenever there is PO chaos at 
some parameter value /i, we expect the periodic orbits to be primarily saddles, 
and if it likewise had virtually uniform PO chaos, then we would expect 
Dim u ([i) = 1. Assuming there are infinitely many periodic orbits, roughly 
half would be regular saddles, with the rest being flip saddles. Furthermore, 
all attracting periodic points are regular. 

Our first goal is to describe the route to (PO) chaos. That is, if at 
/ii there is no chaos, while at fi2 there is virtually uniform chaos, we explain 
what must happen in this interval in order for chaos to arise. 

We believe that generally there is one typical route to chaotic dynamics. 
Namely, there must be infinitely many period-doubling cascades when /i is 
between fii and fi2- (Each of these cascades in turn has infinitely many 
period-doubling bifurcations.) 

Generic maps. Our results are given for generic maps of a parameter. 
Specifically, we say that the map F is generic if all of the bifurcation orbits 
are generic, meaning that each bifurcation orbit is one of the following three 
types: 

1. A standard saddle- node bifurcation. (Where "standard" means the 
form of the bifurcation stated in a standard textbook, such as Robin- 
son [IT].) In particular the orbit has only one eigenvalue A for which 
|A| = 1, namely A = 1. 

2. A standard period-doubling bifurcation. In particular the orbit has 
only one eigenvalue A for which |A| = 1, namely A = — 1. 

3. A standard Hopf bifurcation. In particular the orbit has only one com- 
plex pair of eigenvalues A for which |A| = 1. We require that these 
eigenvalues are not roots of unity; that is, there is no integer k > for 
which A fc = 1. 

These three bifurcations are depicted in Figures [7] and [9] Generic F have 
at most a countable number of bifurcation orbits, so almost every fi has no 
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Figure 9: A regular periodic Hopf bifurcation point is locally con- 
tained in a unique curve of regular periodic orbits. To understand 
generic Hopf bifurcations, it is sufficient to determine the dynamics on its 
(invariant) center manifold, which is locally like R 2 x R. It contains all 
periodic orbits that are near the bifurcation orbit. Near a generic Hopf bi- 
furcation with no eigenvalues which are roots of unity, the local invariant 
set consists of a curve of periodic points (in solid black) and a surround- 
ing paraboloid (at left). Although there are infinitely many periodic orbits 
on the paraboloid in each neighborhood of the Hopf bifurcation point, the 
middle curve of periodic points is disconnected from all periodic points on 
the paraboloid. Specifically, the paraboloid contains infinitely many annuli 
converging to the bifurcation point, each with periodic orbits formed and 
destroyed at saddle-node bifurcations. Typical annuli are depicted here by 
colored regions. The periodic orbits are depicted in blue (saddles) and red 
(nodes). Between every two annuli, there is an invariant circle (in black) 
with irrational rotation number, as made clear by the projection (at right) 
of the paraboloid to the center manifold plane of the phase space. 



bifurcation orbits. See [13] for the details showing that these maps are indeed 
generic in the class of smooth one-parameter families. For systems with 
symmetry such as the forced-damped pendulum, a fourth type of bifurcation 
occurs, such as a pitchfork or symmetry-breaking bifurcation. This adds 
complications, though in fact with extra work our results remain true. 

Our motivation for considering generic maps is given in Proposition [T] in 
the next section, which states that each regular periodic orbit for a generic 
map is locally contained in a unique path of periodic orbits. The connection 
to cascades can be summarized as follows: starting at each regular periodic 
orbit Q for fi e [fii, fx 2 ], there is a local path of regular periodic orbits through 
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Q. Enlarge this path as far as possible. Either the path reaches or /i 2 , or 
there is a cascade. This idea is explained in more detail in the next section. 

3 Onset of chaos implies cascades 

Our first result is Theorem [TJ which demonstrates that the route to virtually 
uniform PO chaos contains infinitely many period-doubling cascades. Theo- 
rem [2] is a restatement of these results in a way that makes the relationship 
between chaos and cascades much more transparent. 

Write J for the closed parameter interval [/ii,/^]- Our main hypotheses 
will be used for a variety of results so we state them here. 

List of Assumptions. 

(A ) Assume F is a generic smooth map; that is, F is infinitely differentiable 
in n and x, and all of its bifurcation orbits are generic. 

(Ai) Assume there is a bounded set M that contains all periodic points 
x) for fi G J. 

(A 2 ) Assume all periodic orbits at n\ and fi2 are hyperbolic. 

(A3) Assume that the number Ai of periodic orbits at fi\ is finite. 

(A4) Assume at /i 2 there is virtually uniform PO chaos. Write A 2 for the 
number of periodic orbits at /i 2 having unstable dimension not equal 
to Dim„(/i 2 ). 

Theorem 1. Assume (A — A4). Then there are infinitely many distinct 
period- doubling cascades between fii and \i2- 

Example 1. Based on numerical studies, a number of maps appear to satisfy 
the conditions of the above theorem. Note that numerical verification involves 
significantly more work than just plotting the attracting sets for each param- 
eter, since we are concerned about both the stable and the unstable behavior 
to determine whether there is chaos. Examples include the time-2^ maps for 
the double-well Duffing (Fig. the triple-well Duffing, the forced-damped 
pendulum (Fig. the Ikeda map (introduced to describe the field of a laser 
cavity), and the pulsed damped rotor map. 
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This formulation gives no idea how the behavior at /i 2 is connected to 
the cascades that must exist. Thus before giving a proof, we reformulate 
the conclusions of the theorem in a more transparent way. Specifically, we 
reformulate so that it is clear how infinitely many cascades in the strip S = 
[111,112] x 9Jt are connected to regular periodic orbits at fi 2 by continuous 
paths in RPO. 

Paths of orbits. We now give a formal definition for paths of regular 
periodic orbits and show that they connect regular periodic orbits at ^2 with 
cascades between and ji 2 . We start with the following theorem, which 
guarantees that a path through each regular periodic orbit is unique. 

Assume hypotheses A — A4, and consider a local path through a regular 
periodic orbit as guaranteed by the above proposition. Specifically, let Y = 
(/i2,xo) be a regular periodic point. Since all such points at /i 2 in the above 
theorem are assumed to be hyperbolic, the orbit can be followed without a 
change of direction for nearby /i < fi 2 . We can define Y(ip) = (/i(ip),x(ip)) 
continuously for ip in an interval [V'ojV'i] as follows. Let Y(if) ) — (^2,^0)- 
Then Y{jp) follows the branch of regular periodic points for decreasing \i 
until reaching Y(ipi), which is either a bifurcation orbit or fj, — fj,i. If it is 
a bifurcation point, which would be a generic bifurcation point, then from 
the above proposition there is a unique branch of regular hyperbolic periodic 
points leading away from Y(ipi), and Y(ip) follows that branch. If Y(ipi) is a 
period-doubling bifurcation point, there are two branches of regular periodic 
points, but both are on the same orbits. Either branch can be followed since 
both represent the same orbits and both are regular. 

Each branch of hyperbolic periodic orbits can be parametrized by \x and 
is easy to find numerically by solving a differential equation. Write A(p) 
for the square matrix D x F p (ii,x(ii)) and b(/i) for the vector D^F p (fi,x(/i)). 
Differentiating 




b(fj.) + A(fj.) 



dx dx 
dji d[i 



dx 
djj, 



(A(v)-Id)- 1 
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Where Id denotes the identity matrix. The matrix inverse above exists 
since (/x, being hyperbolic, cannot have +1 as an eigenvalue. 

This construction suggests a definition. We say Y(ip) is a path in S for 
if) in some interval K if it is continuous and 

(i) Y(ip) is a regular periodic point in S for each ip 6 K; 

(ii) Y does not retrace orbits; that is, Y is never on the same orbit for 
different ip. 

Proposition 1 (Local paths of regular periodic points). For a smooth generic 
map F : R x VJH —> VJt, each regular periodic point is locally contained in a 
path of regular periodic points. The corresponding path of periodic orbits is 
unique. 

This proposition appears in reference [13J. Here we give an idea of the 
proof. The implicit function theorem shows that there is a local curve of 
periodic points through a regular hyperbolic periodic point. Therefore, the 
proof consists of showing that there is also a local curve of orbits through 
each periodic bifurcation orbit. This must be shown for each of the three 
types of generic bifurcations. The curve is not necessarily monotonic: It 
may reverse direction with respect to /x, as depicted in Figure [6j Namely, a 
path Y(ip) = (n(ip),x(ip)) reverses direction at ipo if MVO changes from 
increasing to decreasing or vice versa at that point. For n-dimensional x, the 
study of generic saddle-node and period-doubling bifurcations can be reduced 
to one-dimensional x due to the center manifold theorem. The bifurcation 
orbit has n — 1 eigenvalues outside the unit circle, and these do not affect 
the dynamics on the two-dimensional invariant (/i, x)-plane. Hence here we 
can consider x as being one-dimensional. 

In the case of a saddle-node bifurcation, one branch of orbits is stable 
and one is unstable in R l . Since the other n — 1 eigenvalues do not cross 
the unit circle near the point, it follows that the unstable dimension of the 
orbits in n— dimensions is even for one branch and odd for the other, i.e., the 
unstable dimension k(ip) of Y (ip) changes parity as Y(ip) passes through the 
bifurcation point Y(ip Q ). 

For a period-doubling bifurcation, consider the notation depicted in Fig- 
ure [6] Namely, let (A) denote the period-n orbits with no coexisting period- 
2n orbits, let (B) denote the branch of period-n orbits coexisting at the same 
parameter values with the branch (C) of period-2n orbits. If branch (A) is 
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regular and (B) is flip, then the path Y includes branches (A) and (C), which 
are on different sides of the bifurcation point, and fi(ip) does not change di- 
rection and k{ip) does not change parity. If however (B) is regular and (A) 
is flip, then (i(if)) does change direction and k(ip) does change parity. 

For Hopf bifurcations (see Figure [9]), the path proceeds monotonically 
past the bifurcation point and the unstable dimension changes by ±2, so the 
parity does not change. Hence in all cases, the path changes direction if and 
only if the parity changes. The cases of saddle-node and period-doubling 
bifurcations are straightforward, as depicted in Figure [7j To show that RPO 
is locally a curve at a Hopf bifurcation is trickier than the other two cases, 
since there are in fact infinitely-many periodic orbits near a Hopf bifurcation 
point, but they are not connected to the Hopf bifurcation point by any path 
of periodic points, as shown in Figure [9j The proof follows from the steps 
listed. 

Cascades. We call a regular path Y(i/j) for ip 6 [a, b) a cascade if the path 
contains infinitely many period-doubling bifurcations, and for some period 
p, the periods of the points in the path are precisely p, 2p, 4p, 8p, ■ ■ ■ . As one 
traverses the cascade, the periods need not increase monotonically, but as 
if) — > b, the period of Y(ip) goes to oo. 

We will say a path Y(ip) is maximal if the following additional condition 
holds: 

(iii) Y cannot be extended further to a larger interval, and it cannot be 
redefined to include points of more regular orbits. 



Figure 10 shows an example of a middle portion of a maximal path. 
Figure 11 and [12] show different possibilities for how the maximal path ends. 

Let Orbits(Y) be the set of periodic points on orbits traversed by Y . That 
is, if (/!, x) = Y{jp) for some ■?/>, then ((jl,x) G Orbits(F) and so are the other 
points on the orbit, (/i, F n (fi, x)) for all n. 

Two integers are said to have the same parity if both are odd or both are 



even. Otherwise they have opposite parity. Figures [6] and [10] demonstrate 
why this is a critical idea. Namely, the unstable dimension of a path Y{ip) 
changes parity as ip increases precisely when the path changes directions. 
Thus the parity of the unstable dimension corresponds to the orientation of 
the path. 

Theorem 2. Assume (A — A^). Then the following are true. 
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II 

Figure 10: Paths of regular periodic orbits. We assume here that all 
bifurcation orbits are generic. It follows that a regular periodic point is in a 
path of periodic points. Since each periodic point is in a periodic orbit, we can 
think of the path as being a path of regular periodic orbits. This is a unique 
path of regular orbits (solid). Here we show orbits, not points. Flip orbits 
are depicted by dashed lines, showing that the path is no longer unique if we 
include all periodic orbits. Let Y(ip) be a path of regular (non-flip) periodic 
points parametrized by ip, such that Y(ip) does not pass through the same 
orbit twice. Write |V(^>)] for the orbit that the point Y(ip) is in. It is a 
useful fact that whenever the fi coordinate fi(ip) of Y(ip) changes direction, 
the unstable dimension of Y(ip) changes by one. Hence it changes from odd 
to even or vice versa. The arrows along the path show which way the path 
is traveling as ip increases. 



(Bl) there are infinitely many regular periodic points at ji 2 - 

(B2) For each maximal path Y{jp) = (/i(ip),x(ip)) in S starting from a reg- 
ular periodic point Y = (// 2 ,£o); the set of orbits traversed, denoted 
by OrbitsiY), depends only on the initial orbit containing Yq. That is, 
different initial points on the same orbit yield paths that traverse the 
same set of orbits, so we can write Orbits(Y ) for Orbits(Y). 

(B3) Let Yq = (/j,2,xo) and Y\ = (fi2,xi) be regular periodic points on differ- 
ent orbits. Then Orbits(Yo) and Orbits(Yi) are disjoint. 

(BJf.) Let K denote the unstable dimension of a regular periodic point (fi2,xo)- 
For a maximal path Y(ijj) = (/i(ip),x(ip)) in S starting from (/i2,xo), 
let k(ip) denote the unstable dimension ofY(ip). At each direction- 
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Figure 11: Paths of regular periodic orbits starting from fi2- Assume 
the bifurcations are generic. We use the notation of the above figure. Let 
— (^2,^0) be a regular hyperbolic periodic point. If its maximal path of 
orbits extends to a hyperbolic orbit Y * at /ii, then the path has changed 
directions an even number of times and the unstable dimensions of Yq and 
Y * have the same parity; that is, both are even or both are odd. If the path 
starting at a regular hyperbolic orbit Y\ = (/i2,£i) returns to a point Y* at 
/i 2 , the path changes directions an odd number of times, and so the unstable 
dimensions of Y\ and Y± have opposite parity and in particular the unstable 
dimensions of the two are different. 



reversing bifurcation, k{ip) changes parity; that is it changes from odd 
to even or vice versa. Initially Y(a) = (^2,^0) so initially fi(tp) is 
decreasing and k(a) = K , so initially k(a) + K is even. Hence in 
general fJ,(t(j) is decreasing if K + k(ip) is even and increasing if it is 
odd. 

(B5) Let Y be a maximal path on [a,b] in S and /i(a) = If /J,(b) = H2, 
then fi(ip) is increasing at that point so k(a) + k(b) is even. Hence 
k(a) + k{b) is odd, so k(a) 7^ k{b). 

(B6) There are infinitely many distinct period- doubling cascades between H\ 
and fa- 

(B7) There are at most A = A1+A2 regular periodic orbits at fi2 with unstable 
dimension Dim u (n2) that are not connected to cascades. 

Theorem [2] is a version of Theorem [T] that allows us to show that most 
regular orbits at ^2 are connected to cascade by a path of regular orbits. The 
conclusions of this theorem were shown to hold in [T5] under assumptions 
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Figure 12: An interior path of regular orbits starting from /i2 yields 
cascades. Continuing the assumptions and notations of the above figures, 
assume the hypotheses of Theorem [TJ Only a finite number of paths that 
start at \ii can either reach \i\ or return to Hence there are an infinite 
number of regular hyperbolic orbits at \i2 that yield paths that neither reach 
[L\ nor return to fi2- Such interior paths must contain an infinite number of 
bifurcations. Since generically there are only a finite number of bifurcation 
orbits of each period between ^ and /i 2 , the period of the orbits must tend 
to infinity along the path and it must contain a cascade. Hence there are 
infinitely many cascades between Hi and 

A — v4 4 along with the additional hypothesis that at fi2, there are infinitely 
many periodic regular orbits. However, we are now able to show that this 
last hypothesis is unnecessary because it is automatically true. In fact, under 
assumptions A — A 4 , approximately half of the periodic orbits are regular. 
This proof is quite technical in that it involves topological fixed point index 
theory. To make this treatment more readable, we give the full details in [14]. 

Example 2. Assume phase space is two-dimensional. If at ^ there is a 
transverse homoclinic point then there will be infinitely many saddles and 
exponential orbit growth. Condition A2 is often satisfied if the attractor(s) 
are periodic orbits and there is a Smale horseshoe in the dynamics. 

Example 3. Our numerical studies of the forced-damped pendulum indicate 
that it satisfies the hypotheses of this theorem on [//i,//2] ~ [1-8, 20], [73, 20], [73, 175], 
and [350,175]. Thus there are infinitely-many cascades on each of these in- 
tervals. 

Weakening the uniformity hypothesis. We end this section with the 
conjecture that A4 can be weakened without effecting the conclusions of the 
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theorem. Specifically, let <&(p,K,p) be the fraction of all period-p orbits at 
p that have unstable dimension K. We will say that most periodic orbits 
have unstable dimension K if K, p) — > 1 as p — > oo. We believe that 
for most generic smooth maps and most p, there is some dimension w for 
which this holds. 

Conjecture 1. Assume (Aq — A3) (not including (A4)). Assume that most 
periodic orbits at p 2 have unstable dimension K . Then there are infinitely 
many regular periodic orbits with unstable dimension K , and infinitely many 
of these are connected to distinct period- doubling cascades between pi and p 2 . 

4 Conservation of solitary cascades 

Let J be an interval, and let Y(ip) for ip e J be a maximal path of regular 
periodic orbits containing a cascade. For any point so in the interior of J, 
infinitely many period doublings of the cascade occur either for ip < sq or 
ip > Sq. In the following definition, we distinguish two types of cascades 
based on what happens to the other end of Y(ijj). 

Definition 1. Let F satisfy Aq — A 2 . A cascade is solitary in [pi,p 2 ] if 
the maximal path of regular periodic orbits containing it contains no other 
cascades. In this case, we call the non-cascade end of the maximal path the 
stem of the cascade. A cascade is paired if its maximal path contains two 
cascades. A cascade is bounded in [pi,p 2 ] if its maximal path is contained 
in the interior of the interval and never hits the boundary p = p% or p = p 2 . 
We call a path which stays in the interior of the interval an interior path. 

On a sufficiently large parameter interval, all cascades for quadratic maps 
are solitary, as shown in Figure |4| Figure [5] shows an example of two sets 
of paired cascades. The following theorem shows that the classification of 
solitary versus paired is equivalent to classifying cascades in an interval as 
being purely interior versus hitting the interval boundary. 

Theorem 3 (Solitary and paired cascades). Let F satisfy A — A 2 on [pi, p 2 ]- 
A cascade is paired if and only if the maximal path containing the cascade is 
an interior path. In particular, bounded cascades are always paired. 

The maximal path of a solitary cascade contains a periodic point on the 
interval boundary; i.e. there is a point (p,x) in the maximal path of the 
cascade such that p = pi or p = p 2 . Namely, solitary cascades are never 
bounded. 
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Proof. We give a sketch of the proof. Let J be an interval, and let Y(ip) for 
ip G J be a maximal path of regular periodic points in [//i,/^] containing 
a cascade as ■?/> increases. Fix sq G J, and consider the one-sided maximal 
paths Z±(i(}) C Y(ip) where Z±(ip) = Y(ip) respectively for ip > s ,ip < s , 
ip G J. We have assumed that Z + (ip) contains the cascade. Therefore Z + (ip) 
is an interior path, since otherwise there would not be infinitely many period- 
doubling bifurcations. 

Key Fact: By the methods described in Section [3j a one-sided maximal 
path of regular periodic points is an interior path if and only if it contains a 
cascade. 

Assume the cascade in Y(ip) is paired. Then Z_(ip) C Y(ip) contains a 
cascade, so by the key fact above, Z-{ip) is an interior path, implying Y{ip) 
is as well. 

Conversely, if Y(ip) is an interior path, then Z_(ip) is as well, and by the 
key fact Z_{^) contains a cascade, meaning that the cascades in Z±(ip) are 
paired. □ 



Conservation of solitary cascades. If F satisfies A — A 4 , and A x = A 2 = 
0, then Theorem [I] implies that there is conservation of solitary cascades. 
Specifically, let F be any generic map such that F and F agree at fi\ and H2, 
though they may have completely different behavior inside the interval. Then 
F and F have exactly the same solitary cascade structure. Since the number 
of cascades is infinite, this does not appear to be a meaningful statement. 
However, we can classify solitary cascades by looking at the period of their 
stems, and it is in this sense that their cascades are the same. In Figure [4] 
we label five solitary cascades, one of period three, one of period four, two of 
period five, and one of period six. 

Three rigorous examples of conservation of solitary cascades are encap- 



sulated in the following one-dimensional maps (see Figure 13) 



F(fx, x) = /i — x 2 + g(fx, x) (quadratic), (2) 
-F(/i, x) = fix — x 3 + g(fi, x) (cubic), 
F(fj,, x) = x 4 — 2fix 2 + fi 2 /2 + g(fi, x) (quartic), 

where g is smooth, and for some real positive (3, 

|(?(/i, 0)| < (3 for all fi, and (3) 
\g x (fi,x)\ < (3 for all /i, x. 
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Figure 13: Quadratic, Cubic, and Quartic Maps. The three maps from 
Eqn. |2j are shown for g = and for fi = 24, 15, and 8, respectively. The 
squares depicted are [— 2^//7, 2y / /I] 2 . The maxima and minima are far larger 
than the size of the boxes, resulting in purely uniform chaotic behavior within 
the boxes. This is typical for large \i. As \x increases, these three graphs would 
be stretched vertically. The values of the map F for large \i at the critical 
points are proportional to /i, jj 3 ^ 2 , and /i 2 , respectively, largely unaffected by 
g, which is bounded and has bounded derivative. 



It is straightforward to show that for g = 0, each of these maps has a 
[^1,^2] suc h that there are no regular periodic orbits for F at /^*, and F 
has virtually uniform (PO) chaos at [i* 2 . This leads to the fact that for g 
of the form given, each of the three maps has no regular periodic orbits for 
/ii sufficiently negative, and for \xi sufficiently large has a one-dimensional 
horseshoe map. The conditions on g guarantee that it does not significantly 
affect the periodic orbits when |//| is sufficiently large; in particular it does not 
affect their eigenvalues, so it does not affect the number of period-p regular 
periodic orbits for large Hence all three maps have no regular periodic 
orbits for fi small and have no attracting periodic orbits for /1 large, and for 
sufficiently large all periodic orbits are contained in the set [—2y/Ji, 2 v //I]. 
This leads to the following result: 

Theorem 4 (Conservation of cascades). For each of the functions F in 
Eqn. ^ F is generic for a residual set of g chosen as in Eqn. [3| For each 
of the three examples, the number of stem-period-k solitary cascades for all 
generic F is independent of the choice of g. 

In other words, fix F to be one of the three types in Eqn. [2] Fix any g as 
in Eqn. [3] such that F is generic. Then there exist //£ and /xm such that as 



Cubic Quartic 
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long as fix < hl and /i 2 > //m> F has the same number of solitary cascades 
on [fii, p, 2 ] as occur in the g = case for the interval l/t*,^]- 

As an interpretation of this statement take a set B as large as you like in 
parameter cross phase space, and make g = — (fi — x 2 ) in B, so F = in B. 
Hence the only periodic orbit in B is the fixed point x = 0, so B contains 
no cascades. It might seem that we have annihilated all cascades by this 
process, but the theorem guarantees that every single one of these solitary 
cascades will appear. They are just displaced from their original location, 
moved outside of B. 

The conservation principle works for higher-dimensional maps as well. 
For example, we have shown in [T2] that even large-scale perturbations of the 
two-dimensional Henon map have conservation of solitary cascades. That is, 
writing x = (xi,x 2 ), 

, , / fx + (3x 2 -x\ + g(fj l ,x 1 ,x 2 ) \ , , 

F(/i,Xi,x 2 )= . w v , (Henon) 

V xi + h{^xi ) x 2 ) J 

where /3 is fixed, and the added function (g, h) is smooth and is very small for 
1 1 Xi, X2) 1 1 sufficiently large. See [12J for a precise (and technical) formula- 
tion. Here the added terms cannot destroy the stems which are unchanged 
in the domain where | |(/i, xi, x 2 )\ \ is very large. Each stem must still lead to 
its own solitary cascade. 

If F{ji, xi, . . . , xn) is a set of N coupled quadratic maps such that 

Xi ^ K(n) — x\ + g(xi, . . . , xn), (coupled) 

where g is bounded with bounded first derivatives, and lim/ c _ 5 . = |- 00 Ki(fi) = 
±00, then there is conservation of solitary stem-period-/c cascades. 

5 Off-on-off chaos for paired cascades 

In the previous section, we concentrated on the implications of Theorem [3] on 
solitary cascades, but it also has implications for paired cascades. It describes 
a situation which is quite common in physical systems, namely the case in 
which parameter regions with and without chaos are interspersed, such as 
show in Figs. [2] and [3j Specifically, we give the following definition describing 
the situation where is no chaos at fi\ and /Z3, while at F is chaotic at fi 2 - 
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Definition 2 (Off-on-off chaos). Assume F satisfies A1-A4 on [/ii,// 2 ] and 
on [yU 3 ,/i 2 ], where [i\ < ^2 < 1^3- Then F is said to have off-on-off chaos 
for ^ < iii < fi 3 . 

If F has off-on-off chaos, we can apply Theorem [2] to both [/ii,/i 2 ] and to 
\P2, ^3] and conclude that there are infinitely many cascades in each of the two 
intervals. The following theorem is stronger, in that we also conclude that 
virtually all regular periodic points at fi 2 are contained in paired cascades. 

Theorem 5 (Off-On-Off Chaos). If F has off-on-off chaos on fii < /x 2 < ^3, 

then F has infinitely many (bounded) paired cascades and at most finitely 
many solitary cascades in (//i,/^ 3 ). 

The result follows directly from combining the results of Theorems [2] 
and HI 

Example 4. Our numerical studies indicate that the time-2n maps of the 
forced double-well Duffing (Fig. and forced damped pendulum (Fig. have 
off-on-off chaos, and that it occurs on multiple non- overlapping parameter 
regions. Here 2tt is the forcing period. We cannot prove there are only 
finitely many periodic attractors, though we find very few. These systems 
have no periodic repellers. 



6 Discussion 



Our results in the context of routes to chaos. We now contrast our 
view with the traditional view of many different routes to chaos. People write 
about Routes to Chaos - where chaos does indeed mean there is a chaotic 
attractor. Below we list classes of distinct routes to chaotic attractors. Our 
results say that for generic smooth maps depending on a parameter, there 
is a unique route from no chaos to chaos (which we have defined to mean 
virtually uniform PO chaos) - in the sense of having a chaotic set that need 
not be attracting, for example when there is a transverse homoclinic point. 
Between a parameter value where there is no chaos and a parameter value 
where there is chaos, there must be infinitely many period-doubling cascades. 
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Routes to a chaotic attractor 



1. A chaotic attractor develops where there was no previous horseshoe 
dynamics; this includes what we might call the Feigenbaum cascade 
route. 

2. There is a transient chaos set and a simple non-chaotic attractor (equi- 
librium point or periodic orbit) and that attractor becomes unstable. 
For periodic orbits there are three ways to become unstable. 

3. Like above except the attractor is a torus with quasiperiodic dynamics. 
How many ways can a torus become unstable? We suspect in many 
ways. 

4. Crisis route: as a parameter decreases, a chaotic attractor collides with 
its basin boundary and the chaotic set is no longer attracting. How 
many ways can this happen? Also unknown. 

5. There is a chaotic attractor and a non-chaotic attractor and as a pa- 
rameter is varied, the initial condition being used migrates into the 
basin of the chaotic attractor. 

6. Homoclinic explosions lead to chaotic dynamics. 

However, using the viewpoint described in this paper, there is only one 
route to chaos. 

Relating PO chaos and positive entropy. Consider the following defi- 
nition, related to exponential growth of periodic orbits, alluded to when the 
definition was given: If \fixed{p)\ ~ G p for large p, then we say G is the 
growth factor, or we call logG the periodic orbit entropy. Taking logs of 
both sides, dividing by p, and taking limits, we get 

,. \og\ fix ed(p)\ . , 

lunsup.^— ^ ^ = logG = h. 

p 

We believe that many of the methods used here will lead to a fruitful study 
of periodic orbit entropy. Note that PO chaos occurs exactly when h > 0. If 
there are only finitely many orbits for example, h — 0. If there is at most a 
fixed number k of periodic orbits of each period p, again h — 0, even if there 
may be infinitely many orbits. In many situations, there is a relationship 
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between exponential periodic growth and positive topological entropy, and 
in fact in some cases periodic orbit entropy is equal to topological entropy 
We elaborate below. 

Bowen showed pQ that for Axiom A diffeomorphisms, you can find topo- 
logical entropy exactly by looking at the growth rate for the number of period- 
p orbits as p goes to infinity For Henon-like maps, Wang and Young proved 
that the topological entropy coincides with the exponential growth rate of 
the number of periodic points of period p [15] . Chung and Hirayama (3] show 
that for any surface diffeomorphism with Holder continuous derivative, the 
topological entropy is equal to the exponential growth rate of the number of 
hyperbolic periodic points of saddle type. That is, you can throw away the 
attractors and repellers. 

In terms of interval maps: Misiurewicz and Szlenk [8] proved that if / is a 
continuous and piecewise monotone map of the interval, then the topological 
entropy is bounded above by the exponential growth rate of the number of 
periodic orbits. Building on [5J, Chung [2] showed that if / is C l+a , similar 
results can be obtained if one counts only the set of hyperbolic periodic orbits 
(namely, periodic points x for which | (f p )'(x) l 1 ^ > 1) or the set of transversal 
homoclinic points of a source. 

Non-smooth maps. Continuous maps F which are piecewise smooth but 
not smooth violate our assumptions, but such maps can be thought of as 
the limits of generic smooth maps F n which differ from F only very near 
the discontinuities of F x . In dimension one, F — p + hx + c\x\ gives a rich 
collection of examples obtained by choosing constants b and c in interesting 
ways. See the extensive literature on border collision bifurcations [10]. Such 
maps deserve more discussion than we can give here but we mention two 
examples. 

(1) For the tent map F(p, x) = p — 2\x\, there are no periodic orbits when 
p < 0; there is one periodic orbit, a fixed point at 0, when p = 0, and orbits 
of all periods when p > 0. All orbits are in families of straight line rays that 
bifurcate from and originate at (0,0), the map's only bifurcation periodic 
orbit. In this case, all bifurcations in all the cascades that would exist for 
the approximating generic families F n tend to (0, 0) as n — > oo. 

(2) For the tent map with slope ±p, namely F(p,x) = 1 — p\x\ for 
1 < jj, < 2, a periodic orbit suddenly appears at x = for countably many 
p. For example if p% denotes the smallest parameter with a period-three 
orbit, then the point x = is a period-three point, and only infinitely many 
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periodic orbits whose periods are multiples of three bifurcate from it. All 
of the cascades for those orbits that we would expect for a smooth map are 
collapsed into the single point (/i3,x = 0). 
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